Kondo impurities in nanotubes: the importance of being "in" 
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Transition metal impurities will yield zero bias anomalies in the conductance of well contacted 
metallic carbon nanotubes, but Kondo temperatures and geometry dependences have not been 
anticipated so far. Applying the density functional plus numerical renormalization group approach 
of Lucignano et al. to Co and Fe impurities in (4,4) and (8,8) nanotubes, we discover a huge 
difference of behaviour between outside versus inside adsorption of the impurity. The predicted 
Kondo temperatures and zero bias anomalies, tiny outside the nanotube, turn large and strongly 
radius dependent inside, owing to a change of symmetry of the magnetic orbital. Observation of 
this Kondo effect should open the way to a host of future experiments. 



PACS numbers: 73.63Rt, 73.23.Ad, 73.40.Cg 

Nanotubes provide a rich playground for a variety of 
many body phenomena, in particular quantum trans- 
port between metal leads pQ . Depending on transparency 
of the electrical contact between the nanotube and the 
leads, conduction may range from insulating with strong 
Coulomb blockade for poor contacts [2], to free ballistic 
transport with conductance close to 4e 2 jh when contact 
transmission is close to one [3] H] . Kondo effects in in- 
trinsic nanotubes have been described, either for poor 
contacts [5JE], and/or in connection with superconduct- 
ing leads [7 j, but none of the classic, extrinsic, single- 
atom impurity type. Here we focus on a high transmis- 
sion lead-nanotube-lead contacts, with a single magnetic 
impurity adsorbed inside or outside a metallic nanotube 
segment - an extrinsic case. Conceptually, this should 
constitute a reproducible system, whose conductance can 
be precisely and predictably controlled by standard ex- 
ternal agents such as magnetic field, gate voltage and 
temperature. The ballistic conductance of the four nan- 
otube conduction channels will be altered by Kondo im- 
purity screening, showing up as a zero bias anomaly [8], 
in a way and to an extent which is presently unknown. 
Transport anomalies have long been reported [9] in a tip- 
impurity-deposited nanotube geometry, at tip-impurity- 
metal systems [TOj |TT] ; and at molecular magnetic break 
junctions [12] - systems with very limited atomic and 
structural control. For an atomistically defined system 
like ours, we aim at theoretical predictions that are not 
just generic - as is often the case in Kondo problems - 
but quantitative and precise about Kondo temperatures, 
conductance anomaly widths and lineshapes. For that 
purpose, we need to implement an ab-initio based pro- 
tocol. According to the "DFT+NRG" formulation by 
Lucignano et al. [13 that goal can be achieved by solv- 
ing a custom-built Anderson model whose parameters are 
determined by the first principles derived impurity scat- 



tering phase shifts. As a specific application, we examine 
here Co and Fe impurities adsorbed on the outside sur- 
face of metallic armchair SWNTs. Results are at first 
disappointing, predicting exceedingly small Kondo tem- 
peratures, and tiny conductance anomalies that would 
be hard to observe. When adsorbed inside the nanotube 
however, the same impurities should yield order of mag- 
nitude larger Kondo temperatures, which moreover in- 
crease with decreasing nanotube radius. When inside, in 
fact, the impurity magnetic orbital symmetry switches 
from parallel to perpendicular to the tube axis, causing 
a dramatic increase of hybridization with the carbon tt- 
orbitals, and a corresponding surge of Kondo energy. 

Following Lucignano et al. [13] we first carry out a 
standard spin-polarized density functional theory (DFT) 
electronic structure calculation of the nanotube with 
one impurity; the conduction 7r-electron phase shifts ex- 
tracted from that calculation are used to fix parameters 
of an Anderson model; the model is solved by Numeri- 
cal Renormalization Group (NRG) to obtain Kondo tem- 
peratures; finally, a non equilibrium Green function tech- 
nique (NEGF) yields the conductance near zero bias. We 
choose (4,4) and (8,8) single wall nanotubes (SWNT) 
(Figs. [TJ [2} with a Co or a Fe atom adsorbed at the 
hexagon center in a fully relaxed position alternatively 
outside or inside the tube (details provided in Supple- 
mental Material). 

The impurity projected density of electronic states is 
shown in Fig. [3] 

DOS peaks mark the impurity d-states. Relative to 
the impurity site, states are even (e) or odd (o) under 
reflection across an xy plane (orthogonal to the tube), 
and symmetric (s) or asymmetric (a) under reflection 
across an xz plane slicing the tube lengthwise. Note 
the "magnetic" orbitals, where up and down spins are 
exchange-split below and above the Fermi level respec- 
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FIG. 1. Sketch of the (4,4) (left) and (8,8) (right) SWNTs in 
the yz plane, with an impurity adsorbed in the hollow position 
(either inside or outside). 



IN 



OUT 



5 (^O 

o 




FIG. 2. Spatial distribution of d xy and d xz orbitals on the 
(4,4) SWNT (different colours mean a change of sign of the 
wavefunction).For Co, d xy (well hybridized) is the relevant or- 
bital when inside, and d xz (poorly hybridized) when outside; 
for Fe, all orbitals are relevant. 




FIG. 3. Symmetry resolved PDOS on the impurity atom 
for Co outside (a) and inside (b) the (4,4) SWNT. In the 
first case, orbital d xz is magnetic, while orbital d xy is weakly 
copolarized, and goes in the mixed- valence regime when the 
AIM is solved. In the second case, both orbitals are partly 
polarized; when the AIM is solved, orbital d xy goes to the 
Fermi energy, while orbital d xz is in the mixed- valence regime. 
Results are qualitatively the same on the (8,8) SWNT, where 
however energies differences between d xz and d xy orbitals are 
even smaller. When Fe is considered instead of Co, both 
orbitals are magnetic in all cases. 



tively. In Co/(4,4)-OUT there is a single magnetic or- 
bital d xz with {0,5} symmetry indicating a S=l/2 state 
for Co (3d 9 4s°), and S=l for Fe (3d 8 4s°) on (8,8). Con- 
sider connecting the two ends of a nanotube segment to 
metal leads, and passing a current. If the contacts are 
transparent, ballistic transport along the metallic nan- 
otube will take place through the two bands at Fermi (see 
Fig. ??, Supplemental Material). Left- and right-moving 
electronic states, <pi and <p r give rise in e and o combi- 
nations, (j) e / = (fa ± r )/\/2 to four channels with dis- 
tinct symmetries {e/o, s/a} implying without impurities 
a conductance 4e 2 /h for perfectly transmitting contacts. 
A single impurity will cause each conduction channel to 
scatter onto the impurity orbital(s) of same symmetry, if 
any, giving rise to a scattering phase shift. The (8 x 8) 
unitary ^-matrix is diagonal with eigenvalues e 2 ^ CT in 
the {e/o, s/a} representation, where \i = es,ea,os,oa, 
a =t, -l- The transmission and reflection probabilities 

|^Q!(j| — COS (^eaa ^oacr)? I^aal — sin (S eaa ^00.(7)1 

a = s,a, also relate via the Friedel sum rule Ap aa (E) = 
\ d8oi dE E ^ t° the extra DOS Ap aa induced by the im- 
purity for each symmetry and spin. Phase shifts calcu- 
lated by DFT are then used to determine the parameters 
of an Anderson impurity model for the impurity. For 
each channel, we introduce spin rotation angles defined 
as ^ = 2(6 ^1 — S^). The most general AIM should in- 
clude four scattering channels, i = es,ea,os,oa, and six 



impurity orbitals (one s and five d), a = 1, . . . , 6, hence 
it is of the form 

H = E ( £fc 4a<W + E V *,° ( C L<*a<x + H.c) I 
ika \ a / 

+ ^ ti,kk' ^ik(j C ik' a + Himpi W 
ikk' <j 

where c\ ka creates a spin a electron in channel i with 
momentum k along the tube, d\ a a spin a electron in the 
orbital a of the impurity. Vik, a is the hybridization matrix 
element between conduction and impurity orbitals, which 
is finite only if they share the same symmetry, while U^k' 
describes a local scalar potential felt by the conduction 
electrons because of the translational symmetry breaking 
caused by the impurity. H imp includes all terms that in- 
volve only the impurity orbitals, which, since the orbital 
0(3) symmetry is fully removed by crystal field, can be 
written as 

Hi m p = ^ ^ (^a 
aa 

+ ^2 Uab n a n b + 2 J ab S a • (2) 

a<b 

where n aa = d\ a d aa , n a = Y, a n ^ and J ab < 0, favoring 
a ferromagnetic correlation among the spin densities S a 
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of the different orbit als. The parameters of this Hamil- 
tonian are fixed by requiring them to reproduce in the 
mean field approximation the ab initio DFT shifts, in 
addition to orbital energies [I3j [14] . The AIM hamilto- 
nian ([I]) found in this way is still numerically prohibitive. 
Since our ultimate goal is transport at low temperature 
and small bias, we can neglect orbitals that within DFT 
are either doubly occupied or empty. In this approxima- 
tion there are two active magnetic orbitals for both Co 
and for Fe, either outside or inside the tube. 

The crucial difference between outside and inside arises 
in the nature of the magnetic orbital. For Co outside, d xz 
has os symmetry, lying in the y = plane (see Fig. [5]). 
Its hybridization, 
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the controlling parameter of the Kondo effect, is small. 
When Co is inside on the contrary, the magnetic or- 
bital is d xy , lying in the z = plane, its radial lobes 
(see Fig. [2]) much more hybridized with the tube con- 
duction channels. The orbital switching between out- 
side and inside is due to a reversal of crystal field, and 
to a different hybridization. In the case of Fe instead 
both d xz and d xy are magnetic, and the change in crys- 
tal field does not play any major role. In both Co and 
Fe, the tangential nature of the outside magnetic orbital 
implies no strong dependence of hybridization upon tube 
radius. Conversely, the radial nature of the magnetic or- 
bital gives rise to a large radius dependence when the 
imprity is inside, where it is better "surrounded" espe- 
cially for smaller tube radius. The impurity- nearest car- 
bon coupling V leads to a hybridization width for a (n, n) 
nanotube, T ~ 7rV 2 p/n ocl/n where p is the (radius in- 
dependent) density of states at the Fermi energy. This 
increase of V with inverse radius explains the increased 
coupling of orbital d xzi for example, of Co on tube (4,4) 
(0.087 eV) with respect to(8, 8) (0.058 eV), only partly 
compensated by a slight decrease of V due to the larger 
curvature. As a consequence, Kondo temperatures are 
predicted to decrease exponentially with increasing tube 
radius - so long as higher subbands can be neglected. For 
very large tubes, our single subband model is invalid and 
higher subbands must be taken into account. 

It is worth here discussing, at least qualitatively, the 
limit of zero curvature, graphene. For Co or Fe on 
graphene the orbitals d xz and d xy are degenerate, and 
occupied by three electrons only. This unstable SU(4) 
symmetry can be broken by e.g., spin-orbit [15] . or by a 
Jahn- Teller distortion, both leading to an ordinary SU(2) 
Kondo effect. In either case, Co/nanotube Kondo is basi- 
cally different from Co/graphene. [15 For Fe/nanotube, 
with two electron in two orbitals, the additional Hund's 
rule coupling, of order 1 eV , is larger than the crystal- 
field splitting, hybridization differences, and spin-orbit 
interaction - hence both orbitals should jointly undergo 



Kondo screening. The same conclusion should apply to 
Fe/graphene even though the two orbitals become degen- 
erate, because the Hund's exchange forces the two elec- 
trons to occupy each a different orbital in a spin-triplet 
configuration. 
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FIG. 4. Kondo temperatures as a function of curvature for Co 
(a) and Fe (b); predixted zero bias anomalies are also shown 
for each case. Dots show calculated values from tab. [IJ lines 
are best fits assuming (independently for negative and posi- 
tive curvatures) log(Tx) = a — b/\x\ where x=curvature. The 
shaded area shows the region of low curvature, where our sin- 
gle band approximation breaks down, and additional terms, 
such as spin-orbit coupling, must be included in the Hamilto- 
nian. As a consequence of considering higher subbands, the 
Kondo temperature, that in our model goes exponentially to 
zero, can saturate at a finite value (dotted line in fig. (a), not 
shown in fig. (b); the saturation value is guessed). In fig. (b), 
an arrow shows the range of temperatures where Fe behaves 
like an underscreened impurity. 



The simplified Anderson impurity models just ob- 
tained are solved by standard Numerical Renormaliza- 
tion Group (NRG) [T6] : we adopt a two-band model, 
which takes into account channels os and ea, and orbitals 
d xy and d xz . The approximate Kondo temperatures ob- 
tained in this manner (details in Supplemental Material) 
are given in Table [TJ column 6. We should warn here 
that Kondo temperature are by construction affected by 
a large error, because of their intrinsic exponential de- 
pendence on parameters. With that caveat, we verify, as 
already stated, that Kondo temperatures turn from very 
small when impurities are outside the nanotube, to large 
and radius dependent when inside. The impurity inside 
the nanotube is therefore the geometry which we propose 
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Impurity 


Nanotube 


Position 


Orbital 


r 


T K {K) 


Q 


Co 


(4,4) 


Out 




0.087 


0.6 


-0.03 


Co 


(8,8) 


Out 




0.058 


0.1 


-0.04 


Co 


(8,8) 


In 


z 


0.126 


25 


-0.10 


Co 


(4,4) 


In 




0.380 


600 


-0.11 


Fe 


(4,4) 


Out 


dxz 


0.092 


0.002 


0.01 








dxy 


0.082 


0.3 


-0.02 


Fe 


(8,8) 


Out 


dxz 


0.062 


io- 7 


0.06 








d X y 


0.044 


io- 8 


-0.02 


Fe 


(8,8) 


In 


dxz 


0.081 


10" 5 


0.09 








d X y 


0.134 


0.01 


-0.01 


Fe 


(4,4) 


In 


dxz 


0.126 


0.05 


0.35 








dxy 


0.396 


30 


-0.06 



TABLE I. Kondo orbitals for each system, with broadening T 
and Kondo temperature Tr- 



for experimental verification. 

The zero-bias conductance anomalies and lineshapes 
are calculated by the Keldysh method for non- 
equilibrium Green functions [17 j. Approximating the 
Kondo resonance by a Lorentzian, the conductance for 
a single channel is a Fano resonance g{v) = j^f^^p^^ 
where v = Vb/^k is the dimensionless bias potential, 
q is the shape parameter, and total conductance is the 
sum of even and odd channels. Results are shown in 
Table [T] and Fig. [4] We generally predict in all cases a 
conductance minimum at zero bias (g ~ 0), except in the 
case of Co inside the (4,4) tube, where hybridization is 
so large that the Kondo effect gives way to a frank reso- 
nant level. Finally, the Fe impurity, having two different 
Kondo temperatures can, especially when inside, behave 
as an underscreened impurity in the range of tempera- 
tures for which T K1 «T« T K2 . 

In summary, we applied a DFT + NRG method to 
obtain first principle predictions of the Kondo effect in 
the conductance of well contacted metallic nanotube seg- 
ments hosting a single transition metal impurity atom. 
One first novelty is our claim of ab initio predicting 
power, which has not been common in Kondo problems. 
The main physical surprise is a strong difference between 
outside and inside impurity adsorption; only when in- 
side, the Kondo temperatures are large and increase for 
decreasing radius. The large radius limit, graphene, is 
shown to be an intrinsically different case. Experimen- 
tally, it should be possible to insert an impurity inside 
a metallic nanotube segment, long enough that the fi- 
nite size level discreetness is smaller than Kondo energies, 
and short enough for strong correlations to be negligible. 
The predicted Kondo temperatures and zero bias anoma- 
lies for small radius nanotubes are very substantial and 
should not only be measurable, but should open the way 
to a variety of new observations in a larger variety of 



metallic nanotubes. 
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